Soft-collinear resummation in deeply virtual Compton scattering 
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We derive an all order resummation formula for the deeply virtual Compton scattering (DVCS) am- 
plitude, which takes into account soft gluon exchanges in the non-singlet quark coefficient function. 
We identify the ladder diagrams responsible in light-like gauge for (q 3 log 2 (a; ± £))" contributions. 
The resummed series results in a simple closed expression. 
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PACS numbers: 13.38. Bx, 13.38.Cy, 13.60.Fz 

In the collincar factorization framework the scatter- 
ing amplitude for exclusive processes such as DVCS has 
been shown [H to factorize in specific kinematical regions, 
provided a large scale controls the separation of short 
distance dominated partonic subprocesses and long dis- 
tance hadronic matrix elements, the generalized parton 
distributions (GPDs) 0. The amplitude for the DVCS 
process 
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with a large virtuality q 2 = —Q 2 , factorizes in terms of 
perturbatively calculable coefficient functions C[x, £, a s ) 
and GPDs F(x,£,t), where the scaling variable in the 
generalized Bjorkcn limit is the skewness £ defined as 



Q 2 



(p + p') ■ (q + q') 
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Lessons from the DVCS amplitude at NLO. Let us 
begin with the discussion of the NLO corrections to the 
amplitude for DVCS ([TJ. After proper renormalization, 
the quark contribution to the symmetric part of the fac- 
torized Compton scattering amplitude reads 
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where the quark coefficient function T q reads 0] 
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The first (rcsp. second) terms in Eqs. (JS|) and ^ cor- 
respond to the s— channel (resp. u— channel) class of di- 
agrams. One goes from the s— channel to the u— channel 



by the interchange of the photon attachments. Since 
these two contributions are obtained from one another 
by a simple (x O —x) interchange, we will restrict in the 
following mostly to the discussion of the former class of 
diagrams. 

Eqs. ([5]) and © show that among the corrections 
of 0(a s ) to the coefficient function, the terms of order 
l ° s J±^^ play an important role in the region of small 
(x±£), i.e. in the vicinity of the boundary between the so- 
called ERBL and DGLAP domains, where the evolution 
equations of GPDs take distinct forms. 

Main steps of our analysis. We start our analysis 
with the observation that in the same spirit as for evolu- 
tion equations, the extraction of the soft-collinear sin- 
gularities which dominate the amplitude in the limit 
x — > ±£ is made easier if one uses the light-like gauge 
Pi = py. We argue that in this gauge the amplitude is 
dominated by ladder-like diagrams, see Fig. [TJ In our 
analysis we expand any momentum in the Sudakov basis 
pi, p2, as k = ap\ + /3p2 + k±_ , where p2 is the light- 
cone direction of the two incoming and outgoing partons 



s = Q /2£). In this basis, 



(Pi = V\ = 0, 2 Pl • p 2 

We now restrict our study to the limit x — > +£. The 
dominant kinematics is given by a strong ordering both 
in longitudinal and transverse momenta, according to 



» \X - £ - Pi - fa P n -l 

1*1x1 « l fc i 2 l « • • • « l*i„l < s ~ q 

\ai\ < ■•• < |a„| < 1. 
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This ordering is related to the fact that the dominant 
double logarithmic contribution for each loop arises from 
the region of phase space where both soft and collincar 
singularities manifest themselves. When x — > £ the left 
fermionic line is a hard line, from which the gluons are 
emitted in an eikonal way, with a collinear ordering. For 
the right fermionic line, eikonal approximation is not 
valid, since the dominant momentum flow along p2 is 
from gluon to fermion, nevertheless the collinear approx- 
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FIG. 1: The ladder diagrams which contribute in the light- 
like gauge to the leading a™ ln 2 "(£ — x)/{x — £) terms in the 
perturbative expansion of the DVCS amplitude. The P2 and 
_L momentum components are indicated. The dashed lines 
show the dominant momentum flows along the P2 direction. 

imation can still be applied. 

Finally, the issue related to the ie prescription in 
Eq. ^ is solved by computing the coefficient function 
in the unphysical region £ > 1. After analytical continu- 
ation to the physical region < £ < 1, the final result is 
then obtained through the shift £ — > £ — ie . 

To prepare for further analysis, we define K n , the con- 
tribution of a 7i-loop ladder to the coefficient function, 
as 



K — p l 
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Now we will discuss in details the main points of deriva- 
tion of /„ starting from the simplest case. 

The ladder diagram at order a s . A careful analysis 
shows that the one-loop diagrams of Fig. [2] are not dom- 
inant for x — > £ . Thus, we concentrate on the box di- 
agram, see Fig. [T]with n = 1. Starting from the dom- 
inant part of the numerator of the Born term which is 
fi = —2 jp\, the numerator of the box diagram is 

tr{ M"[>ft+(*-0 M fi[feH*+t) foh v }d» v . (11) 

We first perform the soft gluon approximation : k + (x + 
€)P2 ~ * ( x + &P2- Since the gluon is almost on-shell, 
the dominant contribution in the gluon propagator, d^ v , 
when expressed in terms of gluon polarization vectors, is 
given by d^ v re — J^x e (!\) e (A) • Writing gluon polarization 
vectors through their Sudakov decomposition 
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FIG. 2: The one-loop diagrams which do not contribute to 



the leading 
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terms in the pi ■ A = light-like gauge. 



allows us to define an effective vertex for the gluon and 
outgoing quark through the polarization sum 



J2 e MV - k ±^x) = ( - k l 
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The numerator, (Num)i, is a— independent and reads 
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We emphasize that the term (x + £)//? arises from the 
cikonal emission from the left fermionic line while the ex- 
pression inside the [ • ■ ■ ] accounts for the fact that gluon 
is not soft from the point of view of the right fcrmion (in 
an eikonal treatment it would reduce to (x — £)//?)■ 

Let us now calculate the integral over the gluon mo- 
mentum k, using dimensional rcgularization J d d k — >• 
§ J" da dfi d d ~ 2 k, {k\ = -fc 2 ). The Cauchy integration 
over a, which also determines the limits of the (3 integra- 
tion, gives 
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with the denominators of propagators L\{Ri) on the left 
(right) quark line and S for s-channel quark line 



L\ = a(x + £)s 
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It turns out that the first line in Eq. (|15|) . which corre- 
sponds to cutting through the exchanged gluon, domi- 
nates. Moreover, the relevant region of integration cor- 
responds to small |/3 + x — £|. The j3 and k integrations 
results in our final one-loop expression 

h = --^^\og 2 (a(x-0), (17) 

where we kept only the most singular terms in the x — > £ 
region and have no control of the value of a within our 
approximation. To fix a, we match our approximated 
one-loop result with the full one-loop result This 
amounts to cut the k 2 integral at Q 2 . The ie term is 
included according to the same matching. This leads to 
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A similar expression holds for the u— channel diagram in 
the limit x — > — £ . 

The ladder diagram at order a 2 . Let us examine the 
next order in the perturbative expansion. There are 
many diagrams contributing but it can be shown that 
in the chosen gauge, the double box diagram dominates, 
and moreover that the relevant region of integration is 
where there is a strong ordering between k 2 j_ an d k\±, 
i.e. k 2 ± 3> fcix- Using the same arguments as in the one- 
loop case, one simplifies the numerator for this diagram, 
after taking advantage of strong ki± ordering, as 
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The simplified numerator for the n th order box diagram 
is obtained by generalizing the argument for the two-loop 
case. One gets 



(Num) n = -4s(x + £)' 
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and the denominators of propagators are, for i = 1 ■ • ■ n, 

L i = ai(x + £)s 

s 2 
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Using dimensional regularization and omitting scaleless 
integrals, the integral /„ reads 



The denominators of fcrmionic propagators read (Fig. [T]) 



L\ = ai(x + £)s , R\ 
L\ = a 2 (x + £)s , R\ 
S 2 = -k 2 + Oft + p 2 + x - Os . 



-fc 2 + aiOft + x - £)s . 
-k 2 + a 2 {Pi+ P 2 + x-C)s, 
(20) 



Using Cauchy integration over ot\ and a 2 , we see that 
the dominant contribution comes from the residues cor- 
responding to the two gluon propagators, and reads 

r£—B ri—x-fl-i poo pkl 
h~ / rfft / rfft / d d - 2 k 2 / d d - 2 k x X 2 , (21) 
Jo Jo Jo Jo 

where 
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The integrals over ki and k 2 are performed by assuming 
that the singularities are regularized within dimensional 
regularization. Taking into account the fact that a scale- 
less integral vanishes, we get the following formula 
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The it and the 2£ factor in the argument of the log- 
arithm are fixed by hand, extending the matching used 
at the one-loop level. These prescriptions are not needed 
for the resummation of [a s log 2 (^p — it)] n terms which 
we want to exhibit and are beyond the accuracy of our 
estimate. 

The ladder diagram at order a™ . We now turn to the 
estimation of all log 2n (x — £) terms. We again assume 
the strong ordering (|8l |9]) in k± and a. The distribution 
of the poles generates nested integrals in 
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kf Kk^-(p 1 + --- + p n + x- £)s 

The integrals over k 1 ---k n are performed again in a sim- 
ilar way as in the case of one- and two-loops, resulting 
in 
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where the matching condition introduced in one-loop case 
is extended to n— loops. 

The resummed formula. Based on the results 
Eqs. (fTBl 1231 129|) . one can discuss the resummed formula 
for the complete amplitude. By combining Eq. (|29[) with 
Eq. (dUJ), we get 
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Although the matching condition introduced at one-loop 
still docs not fix the form of the resummed formula 
uniquely, since it requires to go beyond double logarith- 
mic accuracy, it is natural to propose the following two 
forms 
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FIG. 3: ReT q as a function of x, for £ = 0.05 and /iF = 
Q: coefficient function at LO (dot dashed-orange) and NLO 
(solid- black), resummed formulas (|31[) (dashed- red) and (|32[) 
(dotted-blue). For simplicity, e q = 1 . 



with L> 
(C +d) 



2%2V and 



res 
2 



cosh 



D\oe 



2£ 



(32) 



_D2 

2 



(x 



These resummed formulas differ through logarithmic con- 
tributions which are beyond the precision of our study. 

Phenomenological perspectives. In Figs. [3] and [4] we 
show the real and imaginary parts of the resummed coeffi- 
cient function (£ = 0.05 and a s = 0.33, corresponding to 
Q = 2 GeV and Aqcd = -2 GeV) together with the LO 
and NLO result, restricted to the DGLAP region. While 
the NLO corrections have sizable effects with respect to 
LO, Fig. [3] shows that the higher order corrections when 
resummed in our way do not dramatically lead to any 
further change of the coefficient function. This hints to a 
stabilized pcrturbative expansion around x = ±£ . Nev- 
ertheless, Figs. [3] and 0] quantifies the sizable effects from 
the ambiguous choice of resummation formula ([3T1 132]) , 
demanding a higher order analysis. 

Conclusions. We have demonstrated that resumma- 
tion of soft-collinear gluon radiation effects can be per- 
formed in hard exclusive reactions amplitudes. The re- 
sulting formula for coefficient function exhibits a differ- 
ent behaviour than both the Born order and the NLO 
results. A detailed investigation of the effect of resum- 
mation on observables such as the magnitude of Compton 
form factors will be performed elsewhere. A related but 
somewhat different expression should emerge in various 
hard exclusive reactions, such as exclusive meson elec- 
troproduction. The simplest extension of our formula 
will address the crossed case of timelikc Compton scat- 
tering [ij. 

We did not study neither the effects of the running 
of a s nor the case of gluon GPD contributions which, 
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FIG. 4: ImT q as a function of x, for £ = 0.05 and /if = 
Q: NLO coefficient function (solid- black), resummed formulas 



3Tjl (dashed- red) and (|32|) (dotted-blue), with e q 
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although they begin at order a s are expected to be im- 
portant in the small £ regime to be accessible at high 
energies Q . It has been customary to perform resumma- 
tion in terms of Mellin moments |6j . In our case, confor- 
mal Mellin moments may be more appropriate, and we 
will address this question in future work. 
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